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ABSTRACT 

We estimated black hole masses and Eddington ratios (L/LEdd) for a well defined sample of local (z < 0.3) broad line AGN from the 
Hamburg/ESO Survey (HES), based on the Hj8 line and standard recipes assuming virial equilibrium for the broad line region. The 
sample represents the low-redshift AGN population over a wide range of luminosities, from Seyfert 1 galaxies to luminous quasars. 
From the distribution of black hole masses we derived the active black hole mass function (BHMF) and the Eddington ratio distribution 
function (ERDF) in the local universe, exploiting the fact that the HES has a well-defined selection function. While the directly 
determined ERDF turns over around L/L^m ~ 0.1, similar to what has been seen in previous analyses, we argue that this is an artefact 
of the sample selection. We employed a maximum likelihood approach to estimate the intrinsic distribution functions of black hole 
masses and Eddington ratios simultaneously in an unbiased way, taking the sample selection function fully into account. The resulting 
ERDF is well described by a Schechter function, with evidence for a steady increase towards lower Eddington ratios, qualitatively 
similar to what has been found for type 2 AGN from the SDSS. 

Comparing our best-fit active BHMF with the mass function of inactive black holes we obtained an estimate of the fraction of 
active black holes, i.e. an estimate of the AGN duty cycle. The active fraction decreases strongly with increasing black hole mass. A 
comparison with the BHMF at higher redshifts also indicates that, at the high mass end, black holes are now in a less active stage than 
at earlier cosmic epochs. Our results support the notion of anti-hierarchical growth of black holes, and are consistent with a picture 
where the most massive black holes grew at early cosmic times, whereas at present mainly smaller mass black holes accrete at a 
significant rate. 
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1. Introduction 

The observed relations between the black hole mass and the 
properties of the spheroidal galaxy component imply a close 
connection between the growth of supermassive black holes 
(SMBH) and the evolution of their host galaxies. For local galax- 
ies a strong correlation between the mass of th e SMBH and the 
luminosity or mass of the bulge component (iMagorrian et al.l 
[1991 iMarconi & HuntillOOl iHiiring & Rixll2004. as well as 
with the stellar velocity dispersion (e.g. Ferrarese & Merritt! 
2000; Gebhardt et al. 2000; Tremaine et al. 2002; Giiltekin et alj 
20091) have been established. Semi-analytical and numerical 
simulations also show the importance of black hole activity and 
their corresponding SMBH feedback for galaxy evolution (e.g . 



the AGNLF has been use d to gain insight into the growth his 



tory of black holes /e.g. |Soltanl 1 19821: 1 Yu & Tremain^ 120021: 



IMarconi et aT1l2004t lMerlonill2004t IShankar et al.ll2009h . and it 



pi Matteo et al.[ 120051; ISpringel et al.' 2005'; Croton et al."200^ 
Catta neo et alT (2006t [Khalatvan et al. 2008; Booth & Schavg 



20^ . It became clear that the central SMBH of a galaxy and 



especially its growth is an important ingredient for our under- 
standing of galaxy formation and evolution. 

Therefore a complete census of the black hole population 
and its properties is required. Active black holes that will be 
observable as AGN are particularly important to study black 
hole growth. A useful tool to study the AGN population is 
the luminosity function (AGNLF). The observed evolution of 
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became clear that most of the accretion occurs during bright 
QSO phases. But, using the AGNLF alone usually requires some 
additional assumptions, e.g. for the mean accretion rate, and thus 
is affected by uncertainties and degeneracies. Disentangling the 
AGNLF into the underlying distribution functions, namely the 
active black hole mass function (BHMF) and the distribution 
function of Eddington ratios (ERDF), is able to provide addi- 
tional essential constraints on the growth of SMBHs. 

To understand the influence of black hole growth on galaxy 
evolution over cosmic time, first the properties of growing black 
holes in the local universe have to be well understood. Thus, it 
is important to derive black hole masses and accretion rates for 
large, well defined samples of AGN. However, measuring black 
hole masses is much more difficult than measuring luminosities. 
Black hole masses for large samples of AGN can not be mea- 
sured directly, but only estimated, using locally established scal- 
ing relations. 

The best method to estimate Mbh for type 1 AGN is rever- 
beration mapping of t he broad line region dBlandford & McKed 
119821 ; |Peterson|[l993h . Assuming virial equilibrium black hole 
masses can be estimated by Mbh = /^blrAV^/G, where /?blr 
is the size of the broad line region (BLR), AV is the broad line 
width in km/s and / is a scaling factor of order unity, which de- 
pends on the structure, kinematics and orientation of the BLR. 
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Although the physics of the BLR i s still not wel l understood and 
thus a source of uncertainty (e.g. lKrolikll200lh . the validity of 
the virial assumption has been shown by the measurement of 
time lags and line widths for different broad lines in the same 
spectrum dPeterson & Wandell I2OO0I: lOnken & PetersonI l2002t 
lKollatschnvll2003h . 

However, reverberation mapping requires extensive and 
meticulous observations and thus is not appropriate for large 
samples. Fortunately, a scaling relationship has been estab- 
lished betwe en Rm.r and continuum luminosity of the AGN, 
Rblr « ( lKaspi etal.l l2000. 2005[ iBentz et al.l 12006). Thus 
it became possible to estimate Mbh from single-epoch spectra 
for large samples, and has been used extensively in the previ- 
ous years for large AGN samples (e.g. McLure & Dunlot)ll2004t 
' iih 



Vestergaardl2004l Kollmeier et al."2006VNetzer & Trakhtenbrot) 
2007HShen et al.l20 08b: Fine et al. 2008; Gavignaud et al. 2008; 



Trump et al.ll2009l) . 



For the measurement of the line width, different measures are 
commonly used, and it is unclear if one is superior to the others 
for estimating black hole masses. Most commonly used is the 
FWHM, but it has been suggested that the line dispersion cr\ine, 
i.e. the second central moment of the line profile, is a better mea - 
sure of the line width (iPeterson et al.ll2004l: ICoUin et al.ll2006h . 
The line dispersion is more sensitive to the wings of a line and 
less to the core, whereas for the FWHM the opposite is the case. 
An additio nal measure of li ne width used is the inter-percentile 
value (IPV. lFine et ani2008h . 

The application of the virial method to large AGN samples 
allowed the estimation of the active BH MF (McLure & Dunlo n 



2004; Shen et al. 2008 
2008; Kellv et al..,200' 



. jGreene & Ho 2007; Vest ergaard et aki 
Vestergaard & Osmeru 200Ca dataset 



that is perfectly suited to study especially low redshift AGN is 
provided by the Hamburg/ESO Survey (HES). In this paper we 
use a local AGN sample, drawn from the HES, to estimate their 
black hole masses and Eddington ratios, and construct the active 
black hole mass function as well as the distribution function of 
Eddington ratios. 

We first present our data and our treatment of the spectra. 
We estimate black hole masses and Eddington ratios from the 
spectra using the virial method. Next, we determine the active 
BHMF, taking care to account for sample selection effects, in- 
ducing a bias on the BHMF. Thereby, we not only constrain the 
local active BHMF but also put constraints on the intrinsic un- 
derlying distribution function of Eddington ratios. Finally, we 
discuss our results in the context of the local quiescent BHMF 
as well as that of other surveys. 

Thoughout this paper we assume a Hubble constant of Hq - 
70kms ' Mpc ' and cosmological density parameters Q„, = 0.3 
and Qa - 0.7. 



2. The Sample 

The sample of low redshift AGN used in this stud y is drawn from 
the Q SO catalogue of the Hamburg/ESO Survey dWisotzki et alj 
I2000h . For a more detailed descriptio n of the survey and th e sam- 
ple used, see our companion paper dSchulze et al.ll20()9l here- 
after Paper I). Here we only give a short summary. 

The HES is a wide-angle, slitless spectroscopy survey, 
mainly for bright QSOs, carried out in the southern hemisphere, 
utilising photographic objective prism plates. The HES covers 
a formal area of ~9500 deg^ in the sky. After digitisation, slit- 
less spectra in the range 3200 A ^ A ^ 5200 A have been ex- 
tracted from the plates. From these spectra type 1 AGN have 



been identified, based on their peculiar spectral energy distri- 
bution. Follow-up spectroscopy has been carried out to confirm 
their QSO nature. The HES picks up quasars with B ^ 17.5 at 
redshifts of up to z ^ 3.2. 

The Hamburg/ESO Survey yields a well-defined, flux- 
limited sample with a high degree of completeness. The survey 
covers a large area on the sky and the quasar candidate selection 
takes care to ensure that low redshift, low luminosity objects, 
i.e AGN with prominent host galaxies, are not systematically 
missed. As in Paper I, we want to use this wide luminosity range 
at low redshift, which is unique for optical surveys, to study the 
low-redshift AGN population. 

To construct such a local AGN sample we selected all AGN 
from the final HES catalogue (Wisotzki et al., in prep.) that be- 
long to the 'complete sample' and that are located at redshifts 
z < 0.3. The sample contains 329 type 1 AGN. Spectra are avail- 
able for most of the objects from the follow-up observations. For 
five objects, spectra were either missing in our database or they 
were of such poor quality that they were deemed not usable for 
our purposes. Thus our sample is 324/329 » 98.5 % complete 
in terms of spectroscopic coverage. 



3. Measurement of Emission Line Widths 

For the estimation of Mbh for our low redshift AGN sample, 
the broad line width of the H/3, or alternatively the Ho-, emis- 
sion line has to be determined. For the measurement of the line 
widths of the Ha and H/S emission lines we fitted the spectral 
region around these lines by analytic functions, i.e. by a multi- 
component Gaussian model plus continuum. Over this short 
wavelength range we approximated the underlying continuum 
as a straight line. The Ha and H/3 lines are fitted by one, two or, 
if required, by up to three Gaussians, based on visual inspection 
of the fits. Due to the limited resolution of the spectra the narrow 
line component could only be subtracted for a few lines, if a clear 
attribution of one fitting component to a narrow line component 
was possible. Thus a narrow component was only subtracted if 
clearly identified in the fit. Care has been taken to avoid contam- 
ination of the lines by contribution from the [O iii] AA 4959,5007 
A lines and the Fe 11 emission to the H/3 line, as well as from 
[N 11] and [S 11] to the Ha line, by fitting them simultaneously 
with the Balmer lines. For details on the line fitting we refer to 
Paper I. 

We use two different line width measurements, the FWHM 
and the line dispersion for comparison, because there is at the 
moment no consensus which is the most appropriate for the es- 
timation of black hole masses. Both can be easily derived from 
the fit. We then corrected the line widths for the finite resolu- 
tion of the spectrograph. We measured the continuum flux at 
5100 A from the continuum fit to the H/) line region. We cor- 
rected the flux for G alactic extinction, using the dust maps of 
ISc hlegel et"al] ( Il998h . and the extinction law of ICardelli etakl 
(i 19891) and computed the continuum luminosity AL,i(5100 A), 
hereafter Lsioo- 

For the estimation of errors we constructed artificial spec- 
tra for each object, using the fitted model and Gaussian random 
noise, corresponding to the measured S/N. We used 500 realiza- 
tions for each spectrum. We fitted these artificial spectra, fitting 
the line and the continuum and measured the FWHM, the line 
dispersion and the line flux. The error was then simply taken 
as the dispersion between the various realizations. This method 
provides a formal error, taking into account fitting uncertainties 
caused by the noise. Thereby we assume that our multi-Gaussian 




fitting model provides a sufficiently precise model of the true line 
shape. Intrinsic deviations of the line shape from the model will 
increase the error A remaining Fe ii contribution at H/3 might 
also increase the error 

For a subsam ple of 21 AGN also included in the SDSS Data 
Release 5 (DR5; lAdelman-McCarthv et al.l2007h . we compared 
our results to the higher resolution SDSS spectra. We fitted the 
SDSS spectra in the same manner as the HES spectra. The cor- 
relation for cTjing is tight (we found a scatter of 0.07 dex for li/3), 
whereas the scatter in the measurement of the FWHM is signif- 
icantly larger (0.18 dex for HfS). This is at least partially caused 
by the narrow component that can be disentangled better in the 
SDSS spectra. In contrast, o-\i„s is less susceptible to the narrow 
line contribution and thus provides a more precise width mea- 
surement for our sample. We also see evidence for an small un- 
derestimation of the line width compared to the SDSS spectrum, 
especially for narrower lines, with a mean deviation of 0.03 dex. 
This might be caused by the lower resolution of the HES spectra 
compared to the SDSS spectra and therefore the stronger influ- 
ence of the resolution correction. 

A comparison of the continuum lumi nosity, FWHM and 
black hole mass with the quasar sample of IShen et al.l (l2008bh 
is in general agreement with our values for the few objects in 
common. 

The line dispersion is more sensitive to the wings of a line, 
thus to the subtraction of the contaminating lines, i.e. Fe ii and 
[Oiii] for H/3 and [Nii] and [Sii] for Ha. On the other hand, 
the FWHM is more susceptible to the line core, thus to a prope r 
subtraction of the narrow component (see iDennev et al.ll2009h . 
For our data the latter seems to exhibit the larger uncertainty. 
Together with the indic ation that ixiine is a preferable width esti - 
mate over the FWHM dPeterson et al.ll2004t ICoUin et al.ll2006h . 
we decided to use crii„e to estimate black hole masses, and only 
give the results using the FWHM for comparison. 

3.1. Relations between Hp and Ha Line Widths 

We see a well-defined correlation between the line widths (both 
FWHM and crii„e ) of the H/3 and Ha emission lines, as shown 
in Fig. [1] To quantify this relation, we applied a linear regres- 
sion between Ha and HyS in logarithmic units, using the FITEXY 



method dPress et al.ll 19921) . that accounts for errors in both coor- 
din ates. We accounted for intrinsic scatter in the relation follow- 
ing [^emaineetal] (I2OO2I) by increasing the uncertainties until a 
per degree of freedom of unity was obtained. 
We found the following relations for the line widths: 

log(rH„ =0.96 log cTH/;- 0.08, (1) 

log FWHMh„ = 1 . 10 log FWHMh/; -0.17. (2) 

The rms scatter around the best fits are 0.11 dex for crung and 
0. 16 dex for the FWHM respectively. The relations between Ha 
and H/3 line properties using the FITEXY method are shown in 
Fig.ID 

The relation obtained for the FWHM sUght ly de viates from 
the rela tions obtained by iGreene 'feHol (l200l and IShen et al.l 
(I2008ah . showing a stronger deviation from a one-to-one corre- 
lation. This might be due to the lower resolution of our data, thus 
the resolution correction has a stronger effect on the line width. 
This is supported by the slightly larger scatter for our relation. 
The scatter in the relation between the line dispersions is lower 
than between the FWHMs, again favouring crune over FWHM for 
our data. 

The H/3 lines are on average broader than Ha with 
<FWHM(HyS)/FWHM(Ha)) = 1.54 and (cr,i„e(H/3)/crii„e(Ha)) = 
1 .29 respectively. This is lar ger than found in other samples (e.g. 
lOsterbrock & Shuded[T98 2t iGreene & H0II2OOI but in general 
agreement with the physical expectation of an increasing den- 
sity or ionisation parameter of the BLR with decreasing radius. 

4. Results 

4.1. Estimation of Black Hole Masses 

We estimated black hole masses for the AGN using the common 
scaling relationship. The sample of quasars analysed is well in- 
side the ranges in redshift, with z < 0.3, and in luminosity, with 
10"*^ < L5100 < 10"*^ erg s"', over which the scaling relationship 
based on reverberation mapping has been established. So the es- 
timated black hole masses do not suffer from an extrapolation of 
this relationship outside the range for which it is observationally 
tested. 
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Fig. 2. Median composite spectra for 3 luminosity bins, show- 
ing the Can K line region. The upper composite shows the 
high luminosity bin (logLsloo > 44.5), the middle composite 
is for the medium luminosity bin (43.6 < logLsioo < 44.5) 
and the lower composite spectrum shows the low luminosity bin 
(logLsioo < 43.6). The Ca ii K line at 3934 A is indicated as the 
dashed line. 



For the scaling relationship bet ween BLR size and contin- 
uum luminosity we use the values of iBentz et al.l ( l2009h : 



log/?BLR = -21.3 + 0.5191og(L5ioo) , 

with L5100 given in erg s ' and Rblr in light days. 
The black hole mass is thus computed by 



BH 



6.7 -f 



^5100 



0.52 



AV 



1044 ergs- V \km/s 



M, 







(3) 



(4) 



where / is the scale factor and AV the line width used, i.e. the 
FWHM or crjine. We computed black hole masses based on both 
width measurements. We prefer the black hole masses using the 
line dispersion and give the FWHM based black hole masses 
as a reference. For t he line dispersion w e used a scale factor 
/ = 3.85, following ICoUin et al.l (1200 6'). This factor has been 
determined by setting the black hole masses, computed from the 
mean spectrum of the reverberation mapping sample and using 
the line dispersion, to the local Mrh - Ct relation for q uies- 
cent galaxies, similar to the work of lOnken et all (I200I . For 
the FWHM we u sed the common scale factor / - 3/4 (e.g. 
iNetzer et alJll990l) appropriate for a spherical BLR. 

We have not corrected our continuum luminosities L5100 for 
their host galaxy contribution. This might lead to an overesti- 
mation of Mbh for lower luminosity AGN, where the host con- 
tribution becomes significant. To disentangle the host contribu- 
tion high resolution HST imaging is required, which is not avail- 
able for our sample. However, the narrow slit used for the spec- 
troscopy and the AGN selection technique already reduce the 
expected host contribution. Thus, the bias introduced by host 
galaxy contamination is expected to be small, but will lead to 
a systematic effect. 
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Fig. 3. Left panel: Distribution of black hole masses. The black 
solid histogram shows the distribution of black hole masses es- 
timated from cTiine, the red dashed histogram shows Mbh esti- 
mated using the FWHM (and a constant scale factor of / = 
0.75). Right panel: Distribution of Eddington ratios. The his- 
tograms are the same as in the left panel. 



To estimate the degree of galaxy contribution to our AGN 
spectra we used the equivalent width (EW) of the Ca 11 K line 
at 3934 A, because this is the only prominent galaxy absorption 
feature not confused by AGN emission features within our spec- 
tral range. As this feature is only prominent in evolved stellar 
populations, the contribution from a very young stellar popula- 
tion might be neglected. However, low luminosity AGN hosts 
are known to have not particularly blue colours and generally do 
not show extremely youn g stellar populations, but rather indica- 
tions for post starbursts ( Vanden Berk et al]|2006l: iDavies et al.l 
120071) . 



Since the mean signal-to-noise in our spectra is not sufficient 
to accurately measure the Ca 11 EW in individual spectra, we con- 
structed composite spectra for three luminosity bins, depending 
on L5100, shown in Fig.|2l While no Can absorption is detected 
for the highest luminosity composite (log Lg 100 > 44.5), it is 
clearly present in the lower luminosity composites. We measure 
EWs of 0.8 A in the medium luminosity (43.6 < log L5 100 < 
44.5) and of 2.0 A in the low luminosity (logLgioo < 43.6) me- 
dian composite spectrum, respectively. 

To estimate the corresponding galaxy contribution, we used 
model spectra from single stellar population model s with differ- 
ent ages and metallicities jBruzual & Charlotl2003h . The low lu- 
minosity AGN, which will show the strongest host contribution, 
are preferentially spiral galaxies. We modeled them by stellar 
populations with ages between 900 Myr and 2.5 Gyr We added 
various constant AGN contributions to the spectra and measured 
the resulting EWs of the AGN+galaxy spectra. To derive the 
galaxy contribution at 5100 A we assumed a flux ratio of the 
AGN of /5100//3934 = 0.64. An EW of 2.0 A, as measured for 
our low luminosity subsample, corresponds to a host contribu- 
tion to L5100 of 35 - 40 %. This would reduce our black hole 
mass estimate by 0.10 - 0.12 dex. The upper limit we can put 
on the host contribution is ~ 50 %, implying 0.16 dex for the 
Mbh estimation. The medium luminosity subsample shows an 
average host contribution of 15 - 20 %, corresponding to an over- 
estimation of Mbh by 0.04 - 0.05 dex. 

We used these estimates to apply average host corrections to 
the continuum luminosities and thus to the black hole masses. 



Schulze & Wisotzki: The local active black hole mass function 



5 



Although these corrections might be wrong in individual cases, 
for the sample as a whole the host contribution is thereby ac- 
counted for as good as possible for these data. We verified that 
our results are not qualitatively affected by applying or neglect- 
ing this correction. The quantitative change in the results is cer- 
tainly very small. 

The distributions of black hole masses using the FWHM and 
the line dispersion are shown in Fig. [3] The usage of the FWHM 
instead of cr\i„g slightly shifts the distribution to lower values, 
with (logMBe) decreasing from 7.90 to 7.77, and also broad- 
ens the distribution, with the standard deviation changing from 
0.65 dex for crii„e to 0.70 dex for FWHM. 

In the following we will only refer to the black hole masses 
using (Tiine- This width estimate provides a more reliable width 
measurement for our data compared to the FWHM, as discussed 
in Section[3] We have verified that our results are fully consistent 
when using the FWHM instead. 

4.2. Eddington ratios 

To compute the Eddington ratio A - Lboi/i'Edd, which can 
be understood as a normalised accretion rate, we estimated 
the bolometric luminosity from the optical continuum luminos- 
ity Lboi = fiLsim, ap plying a bolometric correction factor of 
fi - 9, as proposed bv lKaspi et alj ( |2000|) . The mean bolomet- 
ric cor rection factor is still somewh at uncertain, ranging from 
7 (e.g. Netzer & TrakhtenbrotI 2OO7I) to va lues around 10 (e.g. 
[Richards et al.ll2006al: iMarconi et al.l 2004), and also seems to 



be de pendent on luminosit v (IMarconi et al.ll2004t iHopkins et al.l 
l2007h. black hole mas s (Kellvetal. 2008) or Eddington ra- 



tio dVasudevan & Fabiaa,2009; .Lusso et al..,20100 . However, as- 
suming a constant value //, at 5100 A is a good approxima- 
tion. The value of = 9 is also in general agreement with the 
value obtained by inte grating over the mean SED presented by 
[Richards et al.l (l2006ah . The Eddington luminosity is given by 
iEdd = 1.3-lO3^M/M0) ergs-i. 

The distribution of Eddington ratios, using the FWHM and 
CTiine, are shown in the right panel of Fig. [3] The mean Eddington 
ratio of this sample is (log A) - -0.92 with standard deviation 
of 0.46 dex using cr\i„e, and (log A) - -0.79 with 0.56 dex devi- 
ation for FWHM. 

This dispersion is higher than that found by other authors in 
higher redshift an d high er luminosity samples dKollmeier et al.l 
[2()06; Shen et al. 2008bl) . Indeed, the shape of the observed dis- 
tribution does depend on the underlying distribution function 
and the selection function of the survey. Thus the observed distri- 
bution of Eddington ratios is affected by the flux limitation of the 
survey and is not a quantity independent of the specific survey. 
The Eddington ratio di stribution will change in mean and disper- 
sion with luminosity ( iBabic et al.l 120071: iHopkins & HernquistI 
I2OO9I) due to this selection effect. Usually it will broaden with 
decreasing typical luminosity. 

This tren d is also cle a rly visi ble in the sample of SDSS AGN 
presented bv lShen et al.l (l2008bh . A redshift dependence is also 
indicated by their data. For their whole sample, covering the 
range 0.1 < z < 4.5, they found a typical dispersion of ~ 0.3 dex, 
similar to the sample of KoUmeier et al. ( 2006) that covers a sim- 
ilar redshift range and includes relatively h igh lum inosity ob- 
jects. Restricting the sample of IShen et aP (l2008bl) to z < 0.3 
gives a deviation of 0.43 dex, similar to our results, but a lower 
mean Eddington ratio of -1.17 in logarithmic units. This trend is 
also present in deeper surveys that cover a wide redshift range. 
In the VVDS a value for the d ispersion of ~ 0.33 dex has been 
found dGavignaud et al.ll2008h . while in the COSMOS survey a 



dispersion of ~ 0.4 dex has been observed (iTrump et al.ll2009h . 
in agreement with our low redshift result. We will discuss this 
issue further in Section|6] 

In Fig. ID we plot black hole mass, Eddington ratio and bolo- 
metric luminosity against each other The first thing we have to 
be aware of when interpreting these plots are the implicit under- 
lying correlations between these quantities. What we effectively 
always show is a combination of continuum luminosity L5100 and 
line width cr\i„f. . Their underlying relation is shown in Fig. |5] 
There is only some week correlation present between Lsioo and 
line width. 

Physically these plots can be understood from the shape of 
the underlying black hole mass function and distribution func- 
tion of Eddington ratios in combination with the selection func- 
tion of the survey, as we will explicitly show in Section |6] The 
Eddington ratio A sp ans the r ange 0.01 < A < 1, similar to 
other optical sfiidies (IWoo & Urrv 2002 iKollmeier et"aD 120061: 
iGreene & Holl2007t IShen et al..,2008bi) . At high values, observa- 
tions have shown that the Eddington rate represents an approx- 
imate upper boundary to the Eddington ratio distribution, im- 
plying a steep decrease of the Eddington ratio distribution func- 
tion toward Super-Eddington values. At low A the sample suffers 
from incompleteness due to the selection effects of the survey. 
This can explain the observed range of Eddington ratios and the 
rough correlation between Mbh and Lboh shown in the left panel 

of Fig. m 

No strong correlation is seen between /land Lboi for this 
low redshift sample. There is a lack of objects in the lower 
right corner of the middle panel of Fig. |4] thus a lack of ob- 
jects with low-A and high luminosity. These objects would have 
Mbh > 2 ■ lO^M© and are rare objects due to the steep de- 
crease of the black hole mass function at the high mass end 
(see Section [STI and |6^ . Thus it is not surprising to see a lack 
of these objects in the sample. The same applies to the lack 
of objects seen in the upper right corner of the right panel of 
Fig- HI These would be objects with relative high Mbh and high 
Eddington ratio. This lack is also caused by the rarity of these 
objects, due to the steep decrease of the black hole mass function 
in combination with the decrease of the Eddington ratio distribu- 
tion function toward the Eddington rate. Therefore, in the local 
universe massive black holes, accreting close to the Eddington 
limit, are exceedingly rare. 

In the right panel of Fig. |4] there is an absence of objects 
in the lower left corner, i.e. objects with low black hole mass 
and low Eddington ratio. These objects are victims of the survey 
selection. They would have low luminosities and therefore only 
the closest would be detectable in a flux limited sample. An ad- 
ditional effect is that the AGN selection in the HES inevitably 
becomes incomplete at Mbj > -19, because the contribution of 
the host galaxy light even to the HES nuclear extraction scheme 
will become substantial, and the object will no longer be distin- 
guished from a normal galaxy, due to a more galaxy like SED or 
due to a no longer detectable broad emission line. Therefore, no 
AGN with Mbj > - 1 8 are detected in the survey, and the range 
-18 > Mb J > -19 is already seriously affected by this survey 
selection effect. Note that lines of equal luminosities in the right 
panel of Fig. |4] are diagonals from the upper left to the lower 
right. This selection effect explains the absence of observed ob- 
jects in this region and results in the apparent anti-correlation 
between Eddi ngton ratio and black hole mass, also seen in other 
samples (e.g. iMcLure & DunlopI l2004t iNetzer & Trakhtenbrol 
2007^. 

In Section |6] we will explicitly show by Monte Carlo sim- 
ulations how the observed distributions arise from an assumed 
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Fig. 4. Left panel: Black hole mass versus bolometric luminosity. Middle panel: Eddington ratio versus bolometric luminosity. Right 
panel: Eddington ratio versus black hole mass. The three lines indicate Eddington ratios of 1 (solid), 0. 1 (dashed) and 0.01 (dotted). 
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Fig. 5. Distribution of the Hy6 line width (criine) with continuum 
luminosity L5100. There is only little correlation seen between 
line width and L5100. 



underlying BHMF and Eddington ratio distribution function un- 
der consideration of the survey selection criteria. 

5. Black hole mass function and Eddington ratio 
distribution function 

5.1. The local active black hole mass function 

The BHMF of quiescent galaxies in the local universe can be es- 
timat ed, based on the rela t ion between Mbh and bulge properties 
(e.g. ISalucci et al.l Il999t IYu & Tremaiii^ 120021: IShankar et al] 
l2004tlMarconi et alJl2004 . Only a small fraction of local black 
holes are currently in an active state, accreting at a significant 
level and appearearing as an AGN. However, AGN do not ac- 
crete at a single value of A, but rather show a wide distri- 
bution of Eddington ratios (e.g. iHeckman et al] 2004 : Yu et al 



2005':'Merlo ni & Heinzl2008l:lHol2009t iKauffmann & HeckmanI 
2009). Therefore it is not obvious what exactly to call an active 
black hole. A pragmatic definition is to use a lower limit for the 
Eddington ratio. A natural choice for such a lower Eddington ra- 



tio for optical type 1 AGN samples would he A ^ 0.01, as this is 
approximately the observed lower value. 

By this definition, our sample suffers from incompleteness at 
low black hole masses, because some low mass and low A AGN 
will be fainter than the flux-limit. The sample is not selected on 
black hole mass or Eddington ratio but on AGN flux. As already 
mentioned, the sample becomes incomplete at Mg, > -19. Thus, 
at low black hole mass only the AGN above this luminosity limit 
will be detected. This introduces a selection effect on the black 
hole mass distribution that needs to be taken into account for the 
determination of the BHMF. In the following, we will refer to 
this selection eff'ect on the black hole mass and the Eddington ra- 
tio distribution as sample censorship, to distinguish it from more 
direct, for example redshift dependent, selection effects on the 
AGN luminosity distribution. 

It is in principle possible to correct for this sample censorship 
by proper use of the respective selection function. If applying the 
usual selection function, which is a function of luminosity, and 
is appropriate for the determination of the luminosity function, 
to the determination of the black hole mass function, incom- 
pleteness is introduced because it has not properly accounte d 
for active black holes below the flux limit dKellv et al.ll2009l) . 
Instead, the selection function has to be derived as a function of 
black hole mass and this selection function has to be applied to 
the construction of the BHMF. However, to do so would require 
knowledge of the, a priori unknown, Eddington ratio distribu- 
tion function. Thus this approach is not feasible without addi- 
tional assumptions. Nevertheless, it can be useful as a consis- 
tency check, as we will show in Section 163] To avoid such ad- 
ditional assumptions, we used a different approach to determine 
the intrinsic underlying active BHMF from our data, taking into 
account the effect of sample censorship. These results are pre- 
sented in Section|6] 

However, in this section we first determine the active BHMF, 
ignoring the effect of sample censorship on the data. We con- 
struct the BHMF using the usual selection function also used for 
the determination of the AGN luminosity function. However, it 
must be kept in mind that in this case we ignore active black 
holes with luminosities below the flux limit of the survey, even 
if their Eddington ratio is high enough to call it active by the 
above definition. Thus, this determined BHMF suffers from in- 
completeness at low mass caused by the sample censorship. 
Nevertheless, this exercise is worthwhile, because it does not 
require any assumptions on the shape of the mass function or 
any information about the Eddington ratio distribution function. 



Schulze & Wisotzki: The local active black hole mass function 
Table 1. Fitting results for the local active black hole mass function, corrected for evolution but not for sample censorship. 
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Schechter 
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DPL (FWHM) 
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Schechter (FWHM) 
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mod. Schechter (FWHM) 


5.15 X 10-" 


6.81 


-0.59 


0.32 


15.33 


1.70 



While the low mass end clearly will be affected by sample cen- 
sorship, the high mass end is already well determined by this 
approach, providing important information on this mass range. 
Also, this uncorrected BHMF can be better compared with pre- 
vious estimates on the BHMF that usually have not properly ac- 
counted for the sample censorship. 

We constructed this active BHMF, not corrected for sam- 
ple censorship, in an equivalent manner as for the determina- 
tion of a luminosity functi on (see Paper I) . We made use of the 
classical 1 / Vmax estimator ( ISchmidtlll968h to construct a binned 
BHMF. The differential BHMF (space density per log Mbh ) is 
thus given by: 



<1)(M.) = 



A log M. ^ 



(5) 



where Vmux is the maximal accessible volume in which the ob- 
ject, with given magnitude could have been found, given the 
flux-limit of the survey and the redshift bin used. The AGN sam- 
ple used has been selected based on UV excess measured in the 
slitless spectra and no selection based on the presence of emis- 
sion lines is applied. Thus, the Vmax values used are equal to the 
ones used for the determination of the AGN luminosity function, 
presented in Paper I. We lack usable spectra for 5 objects, and for 
an additional 7 objects we could not fit H/3 due to poor quality 
of the spectra in this region and/or due to a low H/? contribution. 
Therefore we could not estimate Mbh for 12 objects. We took 
this into account in the survey selection function by multiplying 
the effective area by a factor of 317/329. The exclusion of the 7 
objects without proper H/? measurement may potentially bias our 
results. Therefore we estimated the H/? width from the Ho' mea- 
surement, using Equation [1] and then estimated Mbh for these 
7 objects. Including these objects results in a consistent BHMF. 
Thus, the in- or exclusion of these objects makes no difference. 

To derive the local {z = 0) BHMF we corrected for evolution 
within our narrow redshift bin, < z < 0.3, as described in 
Paper I. We applied a simple pure density evolution model within 
the redshift bin, i.e. p{z) = (1 + z)*" with ko - 5, thus adjusting 
our BHMF to redshift zero. This specific value ensures a result 
of the y/Vmax test consistent with {V/Vmax) = 0.5, as would be 
expected in the case of no evolution. 

The differential active BHMF of the HES is computed for 
bins of A log Mbh = 0.25 dex in the range lO*" < Mbh < 10'' ^ 
The resulting differential local BHMF, not corrected for sample 
censorship, is shown in Fig. |6l 

We used the following functional forms to fit the BHMF. A 
double power law, given by: 



<^(M.) 



*7M. 



(M./M,)-" + (M./M,)-P 



(6) 
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Fig. 6. The differential active black hole mass function for z = 0, 
not corrected for sample censorship. Filled black symbols show 
the BHMF using the line dispersion to estimate the black hole 
mass. The dashed line shows the double power law fit to the 
BHMF, the dotted line gives the Schechter function fit and the 
dashed dotted line represents the fit using a modified Schechter 
function. 



where M» is a characteristic break mass, (p* the nor malisation 
and a and /3 are the two slopes. A ISchechterl (Il976h function, 
given by: 



'M. 

M. \mI 



exp 



M. 



(7) 
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is also used. 

We additionally used a functional form, motivated by the 
quiescent BHMF. The quiescent BHMF is given as a convolu- 
tion of a Schechter function with a Ga ussian and can be param- 
eterised by the follo wing function (e.g. lAller & Richstonell2002l 
IShankar et ani200l : 



(8) 



This basically corresponds to an ad-hoc modification of the 
Schechter function with an extra parameter fi. A value /3 > I 
corresponds to a decrease stronger than exponential and < I 
corresponds to a milder than exponential decrease. Fot/3 = 1 this 
function turns into the usual Schechter function. In the following 
we refer to it as the modified Schechter function. 

These BHMFs are connected to the expression in logarithmic 
units by <[>(M.) = (M./ logm e)0(M.). The resulting fitting pa- 
rameters of these three functions to our binned BHMF are listed 
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Table 2. Binned black hole mass function, not corrected for sam- 
ple censorship. N gives the number of objects in each bin, log (p 
and A log (p gives the space density per unit logarithmic black 
hole mass in solar masses and its Icr error respectively. 
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in Table [T] All give acceptable fits, while the Schechter func- 
tion performs poorly at the highest black hole masses. However, 
the BHMF is less well constrained at high Mbh due to the small 
number of objects in these bins. 

The shape of the BHMF is described by a steep decrease of 
the space density towards higher Mbh with yS -3 in the double 
power law, and a significant flattening at Mbh ~ lO^M© toward 
lower Mbh ■ The high mass regime is not affected by the already 
mentioned sample censorship, while the low mass flattening is 
partially caused by the systematic underrepresentation of low X 
objects at low mass. 

5.2. The local Eddington ratio distribution function 

Given the estimates of the Eddington ratio A for our sample, 
we can analogously determine the local Eddington ratio distri- 
bution function (ERDF) for the HES, equivalent to the BHMF. 
This determination also does not take into account the effect 
of sample censorship. We computed the local ERDF in bins of 
Alog/i = 0.25 dex in the range -2.25 < log.! < 0.25. The 
resulting differential local ERDF is shown in Fig.Q 

The uncorrected AGN space density declines at high as well 
as at low A , showing a peak around log A ^-1.0. We fitted the 
ERDF by a Schechter function, neglecting the lowest A point. 
The resulting best fit values are 0.i * = 6.66 X 10-'' Mpc-\ 
log A, = - 1 .01 and a a = -0.05 with a value ofx^ per degree of 
freedom of 1.9. 

However, also this ERDF is strongly affected by sample cen- 
sorship. While at the highest Eddington ratios {log A > -1) 
the majority of AGN will be detected by the survey, at low 
Eddington ratio (log A < - 1 ) a significant number of objects will 
have a too low luminosity to be detected. Therefore the space 
density at low A will be underestimated by the derived ERDF. 
In the next section we will reconstruct the intrinsic underlying 
ERDF as well as the intrinsic BHMF. 




og X 



Fig. 7. The differential Eddington ratio distribution function for 
z = 0, not corrected for sample censorship. The dashed line 
shows the best Schechter function fit. 

6. Reconstruction of the intrinsic BHIVIF and ERDF 

6.1. Method 

As already noted, the BHMF presented so far is basically lumi- 
nosity limited and thus incomplete at low mass in terms of an 
accretion rate limited active BHMF. We now want to constrain 
the intrinsic active BHMF by our observations, correcting for 
this sample censorship. We use log A = -2 as the lower limit for 
the Eddington ratio; for log A = -2 we call black holes 'active'. 

The selection function of the survey is a function of luminos- 
ity, and thus of the product of Mbh and A . Therefore, the recon- 
struction of the active BHMF also requires the knowledge of the 
ERDF. Both distribution functions cannot be determined inde- 
pendently from each other In Section 1631 as a consistency test, 
we will determine the active BHMF assuming a specific ERDF. 
But without such an assumption both distribution functions have 
to be determined at the same time. This is the approach we will 
follow in this section. 

Knowing both distribution functions, the AGN luminosity 
function is directly given as a convolution of the two: 



0(L) 



Jm.„ 



Pa iA)^.{M.) dlogM. 



(9) 



where we adopt log M, „ 
malised ERDF, thus: 



Pa W 



J (A)dlogA 



6. With Pi we define the nor- 



(10) 



We determined the BHMF and ERD F together, performin g 
a maximum likelihood fit to the data (e.g. [Marshall et aLlll983h . 
We consider the joint Poisson probability distribution of black 
hole mass and Eddington ratio. We minimise the function S = 
-2 In X, with X being the likelihood of finding the observed data, 
given the respective model. Thus, S is given by: 

N „„ 

5 = -2^1np(M.,„i,) + 2 11 /?(M.,i)dlog.ldlogM. .(II) 

The sum is over the observed objects and the integral is equal 
to the expected number of objects, given the assumed BHMF 
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and ERDF. The probability distribution p(M, , /I )d log /I d log M, 
gives the probability of finding an AGN with black hole mass be- 
tween log M, and log M, +d log M, and Eddington ratio between 
log A and log A +d log A in an observed sample. The total num- 
ber of objects is then given by integration of p{M,,A ) over 
M, and A . 

We will now briefly motivate the used probability distribu- 
tion p{M,,A) for our sample. The observed number of objects 
in a sample is given by: 



ff 



dV 

Qeff (m, z) OiQog L,z) — d log L dz , 
dz 



(12) 



where <I>£(log L) is the AGN luminosity function and Qgff is the 
effective survey area as a function of apparent magnitude and 
redshift, thus the selection function for our flux-limited survey. 
This can be understood as a selection function depending on z, 
M, and A, thus Qeff(OT,z) - Qeif(L, z) = Q.eff{M,,A,z)- For the 
Hamburg/ESO Survey the selection function within our covered 
redshift range is almost independent of redshift. Thus, as dis- 
cussed in Paper 1, we can marginalise o ver redshift. For detai ls 
on the selection function of the HES see lWisotzki et al.l (l2000h . 

Apart from the flux limit, our sample is incomplete at the 
lowest luminosities Mbj > -19. For low luminosity AGN the 
host galaxy contribution becomes an important factor and the 
objects might no longer be classified as an AGN, due to the SED 
being dominated by starlight. As shown in Paper 1, the sample 
is highly complete brighter than Mb, ~ -19. Thus, we adopted 
a luminosity limit of Mbj < -19 in the selection function, Qeff. 
We also restricted the observed sample to this lower luminosity 
for the comparison of the sample properties. 

The AGN luminosity function <t)i(logL) is related to the 
BHMF and the ERDF via Equation |9] For the redshift evolu- 
tion, we assumed the simple pure density evolution model of 
Section 15.11 In this case the black hole mass function is sep- 
arable into a function of M, and a function of z, <t>(M,,z) - 
(D(M.)p(z) = 0(M.)(1 + zf", with ko = 5. 

The expected number of objects for a given survey and an 
assumed BHMF and ERDF is then given by: 



fieff Pa (A)(D.(M.) (l+z)*''^dlog AdlogM.dz .(13) 
dz 



Thus, the bivariate probability distribution of black hole mass 
and Eddington ratio is given by: 



p{M.,A) 



^,sPAm.(M.){\+zf''^dz. 

dz 



(14) 



Given this bivariate distribution for an assumed BHMF and 
ERDF, we minimise the likelihood function S (Equationpn i us- 
ing a downhill simplex algorithm (iNelder & Meadlll965l) . As a 
lower limit for the fitting we employed a black hole mass of 
A/min = IO^Mq and an Eddington ratio of = 10-2. The 
HES sample was restricted to these limits accordingly. 

For the BHMF we assumed three different models. Firstly we 
used a double power law with the high mass slope fixed to the 
value ySee = -3.01, determined from the uncorrected BHMF in 
Section lSTI This lowers the required number of free parameters 
and is justified, because the high mass region in the uncorrected 
BHMF is only weakly affected by incompleteness. Secondly we 
also used a double power law, but leaving the high mass slope as 
a free parameter, to be determined in the fit. As third model we 
used the function given by Equation[8] thus a modified Schechter 
function. The starting values for the minimisation algorithm are 
taken from the fit to the uncorrected BHMF. 



We decided to model the ERDF by a Schechter function, 
corresponding to an exponential cutoff close to the Eddington 
limit and a wide power law-like distribution at low Eddington 
ratio. This parameterisation differs from the often assumed log- 
normal distribution. However, a log-normal distribution is only 
motivated by the observed distribution, not accounting for any 
selection effects. Also, a log-normal distribution enforces a max- 
imum and a turnover at low A . A Schechter function is more 
flexible, allowing for a turnover at low values, but not enforc- 
ing it. In particular, it allows an increase of the space den- 
sity at low A . This shape would be consistent with observa- 
tions oft ype2AGN (Yu et al. 2005: Hopki ns & Hernquistl2009t 
Kauflj'mann & Heckman 2009 ), with estimates for the total AGN 
population (jMerloni & Heinzl l2008l) as well as with model ex- 
pectatio ns of AGN lightcur ves from self-regulated black hole 
growth (lYu & Lul2 008: Hop kins & He rnauist 2009). Aside from 
the Schechter function parameterisation of the ERDF, we addi- 
tionally tested a log-normal ERDF as functional form. Together 
with the Schechter function it covers a wide range of possible 
parameterisations for the ERDF. 

From our data we are not able to constrain a dependence of 
the ERDF on Mbh , so we assumed the ERDF to be independent 
of Mbh , already implicitly assumed in Equation|9] The normal- 
isation of the ERDF is fixed by the condition that the BHMF 
and ERDF have to predict the same space density of AGN. This 
leaves two free parameters for the ERDF, the break A , and the 
low-zl slope for the Schechter function, or the mean A , and 
the width cr^ for the log-normal distribution. However, these two 
parameters in both cases are not independent from each other, 
because the data by construction needs to be consistent with the 
observed luminosity function (LF). Thus, for a given BHMF and 
a fixed value for A », is given by the condition that the LF de- 
rived from the BHMF and ERDF by Equation|9]has to be consis- 
tent with the observed LF. Our approach automatically ensures 
the consistency of the BHMF and the ERDF with the observed 
LF 

To assess the goodness of fit for the individual models we 
used two different methods. This is required because the maxi- 
mum likelihood method does not provide its own assessment of 
the goodness of fit. First, we used a two-dimensional K-S test 
(iFasano & Franceschinil IT987I) on the unbinned data. Second, 
we employed a test, binning the data in bins of 0.5 dex in 
Mbh and A respectively. The results are given with the best fit 
parameters in Table [3] 

6.2. Results 

The first model consists of a double power law BHMF, with the 
high mass slope fixed to /3 - -3.11, and a Schechter function 
ERDF. The best fit distribution functions are shown as black 
dashed line in Fig.[8]and their fit parameters are given in Table[3] 
The BHMF shows a steep high mass slope with q-bh ~ -2 and 
the break is consistent with the uncorrected BHMF. The ERDF is 
increasing towards low A down to the applied limit of A = 0.01. 

This function provides a good fit to the high mass end of the 
uncorrected BHMF, which is only little affected by sample cen- 
sorship. At the low mass end the uncorrected BHMF strongly un- 
derpredicts the active black hole space density, compared to the 
reconstructed underlying active BHMF. This also holds true for 
all other applied functional forms for the BHMF and the ERDF. 
The same also applies to the uncorrected ERDF. The uncorrected 
ERDF is strongly biased and underestimates the BH space den- 
sity. The best fit to the uncorrected BHMF and to the ERDF is 
clearly rejected by the maximum likelihood approach with high 
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Table 3. Fitting results for the active black hole mass function and the Eddington ratio distribution function. The first column 
indicates the function used for the BHMF. 'DPL' is for a double power law, with (3 indicating the fixing of the high mass slope and 
'mS' is for a modified Schechter function. The second column indicates the ERDF. 'S' is for a Schechter function and 'In' stands 
for a log-normal distribution. 
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Fig. 8. Results for the reconstructed BHMF and ERDF. The left panel gives the BHMF and the right panel the ERDF respectively. 
The black points show the binned uncorrected distribution function, with filled circles representing bins that do not suffer signifi- 
cantly from sample censorship and open circles represent bins, biased by sample censorship. They are shown for comparison with 
the reconstructed BHMF and ERDF. The black dashed line shows a double power law BHMF with fixed high mass slope p = -3.11 
and Schechter ERDF, the blue dashed dotted line is for a free double power law BHMF and Schechter ERDF, the red solid line 
represents a modified Schechter BHMF and Schechter ERDF and the green dashed line is for a modified Schechter BHMF and log- 
normal ERDF. The gray areas show the Icr confidence regions of both distribution functions, taking into account all 4 parametric 
models. 



confidence. They are not able to produce the observed distribu- 
tions of /V/bh and A and are not consistent with the AGN LF. 
This clearly shows that the usual approach used to construct an 
uncorrected BHMF and ERDF is strongly biased. 



We briefly want to illustrate how the maximum likelihood 
approach is able to reject certain models for the BHMF and 
ERDF and favour others. To compute the expected distributions 
within a grid of free parameters, we restricted the number of 
parameters to two. We fixed the break and normalisation of the 
BHMF. Thus, with the high mass slope already fixed, the only 
free parameter for the BHMF is the low mass slope (Xbh- For the 
ERDF there are two free parameters, the break and the low-/l 
slope of the Schechter function. However, one of these is fixed 
by the constraint to recover the observed AGN LF. We took a a 
as a free parameter and determined the break by ?lx' minimisa- 
tion of the LF computed via Equation|9]to the observed LF. The 



normalised observed distribution of log M, and log A are given 
by: 



ogM.)=l| 



p{M.,A)d\ogA 



p{\ogA) 



p(M.,A)d\ogM. 



(15) 



(16) 



For illustration, in Fig.|9]we compare these expected distri- 
butions with the observed ones within a grid of free parameter 
ffBH and a,x . For a too steep BHMF the number of low mass ob- 
jects is larger than observed, while for a too flat BHMF the num- 
ber of low mass objects is lower than observed. A steep ERDF 
corresponds to a break of the ERDF close to the Eddington limit, 
thus more objects above the Eddington limit and less at low A 
are predicted, compared to the observations. For a too flat ERDF, 
the break needs to be at a low value of A and thus not enough 
objects close to the Eddington limit are predicted. 

We also carried out Monte Carlo simulations for a grid of 
free parameters q-bh and , using the same assumptions as 
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log A4h log >' log A4h log A log Mbh log A 



Fig. 9. Comparison of the expected distribution of Mbh and A (solid lines) with their distribution in the HES sample (red histogram) 
for 9 different combinations of the two free parameter obh and a,i . The central panel is close to our best fit solution. 



above, as well as for the best fit model of the maximum likeli- 
hood estimation. Here we proceeded as follows: First each AGN 
gets assigned a redshift, then its black hole mass is drawn from 
the assumed BHMF, and finally an Eddington ratio is drawn 
from the ERDF. From these values absolute and apparent Bj 
magnitudes are computed, applying a bolometric correction. By 
means of the apparent magnitude By it is decided if the object is 
selected by the survey or not, taking into account the flux-limit. 

We ran Monte Carlo simulations for a wide range of ubh and 
a,i and found results consistent with what we discussed above 
and what is shown in Fig. |9] The Monte Carlo simulations are 
clearly able to discriminate between models that are consistent 
with the data and those that are not. The best matching solutions 
of the Monte Carlo simulations are consistent with the best fit 
from the maximum likelihood method, although 'best matching' 
is not as well defined in this case. 

In Fig. [To] we show the mean of 10 Monte Carlo realiza- 
tions of this best fit model. We show the observed distributions 
for the sample for this model as well as the uncorrected BHMF 
and ERDF, as well as the Mgj -LF and bolometric LF that would 
be determined from an 'observed' sample. To construct such 
an 'observed' sample we again limited the simulated sample 
to Mbj < -19. In the middle panels of Fig. (TO] we then com- 
pare these expected distribution functions with the uncorrected 
BHMF and ERDF determined with the same restriction applied 
(shown as open red symbols). The distributions as well as the 
constructed distribution functions are consistent with the ob- 
served distributions and distribution functions. For models that 
are found to be not consistent with the observations based on the 
maximum likelihood approach, the Monte Carlo samples also 
provide a poor match to the observed distributions and distribu- 
tion functions, and thus can also be rejected based on the Monte 
Carlo simulations. 



These Monte Carlo simulations show that the observed dis- 
tribution of objects between LboI, ^bh and A , as shown in the 
plots of Fig. m are well understood by the underlying BHMF 
and ERDF and the selection function of the HES. These results 
do not qualitatively change using a different functional form for 
the BHMF or ERDF 

As a second model we again used a double power law, but in- 
cluded the high mass slope as an additional free parameter to be 
determined in the maximum likelihood fit. The result is shown as 
blue dashed dotted lines in Fig. [8] The BHMF is highly consis- 
tent with the previous result, with a mild steepening of the high 
mass slope when this parameter is allowed to change in the fit. 

Third, we also used the function given by Equation[8] thus a 
modified Schechter function. The best fit result is consistent with 
the double power law fit over most of the mass range and only 
decreases stronger at the high mass end. All three models are 
good representations of the observed data and therefore span the 
range of acceptable distribution functions. Formally, the modi- 
fied Schechter function has the lowest value of S and the highest 
probability both in the KS-test as well as in the ;^'^-test and we 
will use it in the following as our reference model. 

Apart from the Schechter function for the ERDF, we addi- 
tionally tested a log-normal distribution. This distribution func- 
tion also provides a good representation of the data. In Table [3] 
and Fig. [8] we give a model with a log-normal distribution for 
the ERDF and a modified Schechter function for the BHMF. 
While the BHMF is nearly unchanged, the ERDF deviates from 
the Schechter ERDF at the highest and lowest values, while be- 
ing consistent over a wide range in between. When enforcing a 
turnover in the ERDF, using a log-normal distribution, the data 
are consistent with such a turnover at low A (log A x -1.8). 
However, there is no evidence for a turnover at higher A , where 
the maximum in the observed Eddington ratio distribution is 
present (log /I a -1). 



12 



I I I I I I I I I I I I I I I I I I I i| I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




Fig. 10. Results of 10 Monte Carlo realizations for the best fit model with an assumed double power law with fixed high mass slope 
for the BHMF and a Schechter function parameterisation of the ERDF. Upper panels; Comparison of the distributions of Mbh , ■I 
and Lboi between simulated sample (black, solid histogram) and observed sample (red, dashed histogram). Middle panels: Same 
as Fig. m but for one simulated sample. Lower panels: uncorrected BHMF and ERDF, Mb, luminosity function and bolometric 
luminosity function. The results for the simulated sample are shown as filled black points. The solid black line shows the true input 
function and the dotted lines show the best fit to the uncorrected distribution functions of the observed sample. The open red circles 
in the BHMF and ERDF plot indicate the individual bins for the observed uncorrected BHMF and ERDF with a restriction of 
Mb, < -19 applied. 
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The log-normal fit indicates rather a flattening of the ERDF 
at the low-A end then a real turnover, because it is cut off be- 
fore the turnover, enforced by a log-normal fit, becomes evi- 
dent. However, the low-A regime is dominated by high mass 
black holes. If there is a mass dependence in the ERDF and the 
ERDF flattens towards high M^h , this would be most promi- 
nent at low A . Such a flatten ing would also be consistent with 
[Ho pkins & Hernquisll ( l2009h . who found evidence for a mass 
dependence in the ERDF of type 2 AGN, with a flatter low A 
slope at high Mbh • 

We take into account the log-normal ERDF in the uncertainty 
range of the determination of the BHMF and ERDF. Formally it 
has a higher probability in the applied statistical tests than the 
Schechter function. However, as mentioned, the main deviation 
compared to the Schechter function is above the Eddington limit 
and close to the lower limit atA- 0.01. The number statistics in 
this regions are low and thus a clear discrimination between the 
two models is not possible. Thus, the Schechter function and log- 
normal distributions indicate the range of acceptable ERDFs. 

We derived uncertainties in the BHMF and ERDF by ran- 
domly modifying the best fit parameters for each model and 
computing the likelihood function S . Using AS - S - Smin for 
each random realization, we c onverted into confidence val- 
ues as suming ax^ distribution ( iLampton et al.|[T976tlPress et al.l 
Il992h . For all models within a certain confidence interval the 
BHMF and ERDF is computed and these functions then span 
the confidence range of the two distribution functions. The total 
uncertainty of the BHMF or ERDF is then the sum of the con- 
fidence ranges of the individual models. In Fig.|8] we show this 
sum of the Icr confidence values for the two distribution func- 
tions as the gray shaded areas. 

So far, we assumed the estimated black hole mass to be equal 
to the true black hole mass. However, this is probably an over- 
simplification. It is known that there is a considerable uncer- 
tainty in Mb h estimates using the virial me thod, probably of or- 
der 0.4 dex dVestergaard & Petersonll2006l) . Accounting for this 
uncertainty might change the reconstructed BHMF and ERDF 
in shape as well as in normalisation. We will investigate this im- 
portant point in detail in future work. 

6.3. BHMF and ERDF from the 1/Vmax method 

As mentioned in Section 15.11 there is also a different approach 
to determine the intrinsic BHMF and ERDF, namely using the 
1 /Vmax method, but directly accounting for the selection effects 
in terms of black hole mass or Eddington ratio completeness im- 
posed on the sample by the AGN luminosity selection. In this 
case, the BHMF and ERDF cannot be determined jointly. When 
using the l/Vmax method the selection effects need to be ac- 
counted for in the determination of the accessible volume of the 
individual AGN, given by: 



ileft —dz ; 



(17) 



where Ogfi is the effective survey area as a function of appar- 
ent magnitude, thus the selection function for our flux-limited 
survey, depending on z, Mbh and A , is Qefj (m) - il^siL, z) - 
^eff{M,,A,z). While for the determination of the luminosity 
function the proper selection function to compute Vmax is given 
by Q.eff{L, z), using it for the determination of the BHMF (as we 
did in Section |5?T] ) will lead to the presence of sample selection 
effects, and thus to the observed underestimation of the space 
density at low mass. 



This bias on the determined BHMF can be avoided by using 
a black hole mass selection function, given by: 



Oeff(M, 



P,,(^)QefF(M.,^,Z)dl0g^ , 



(18) 



where (A ) is the normalised ERDF, given by Equation \TU\ 
However, this approach requires knowledge of the ERDF as 
prior information, which is not present a priori. 

Likewise, the ERDF can be derived in an unbiased way by 
using the Eddington ratio selection function for the survey, given 
by: 



PBH(M.)Qeff(M.,^,Z)dlogMBH , 



(19) 



where Pbh(A^bh) is the normalised BHMF, similar to P^iA), 
defined by: 



Pbh(M.) 



0.(M.) 



/(l).(M.)dlogM. 



(20) 



This reqires knowledge of the BHMF, which is also unknown 
beforehand. Thus this approach is usually not feasible for the 
determination of the intrinsic BHMF and ERDF directly from 
the data. 

However, this approach has the advantage that no prior as- 
sumptions on the shape on the ERDF are required for their de- 
termination, once we fixed the assumed BHMF. The same is 
equally true for the determination of the BHMF. The only nec- 
essary information beforehand is on the shape of the ERDF. The 
problem is that one distribution function needs to be known to 
determine the other one. 

Nevertheless, first we can use it as a consistency test, con- 
structing the BHMF from the constraints on the ERDF from 
Section 16.21 and vice versa. The resulting binned BHMFs and 
ERDFs using the 4 best fit models are shown as filled sym- 
bols in Fig. [TT] together with the best fit solutions to the ac- 
tive BHMF and the ERDF, as determined in Section These 
binned BHMFs as well as the binned ERDFs for all 4 models are 
fully consistent with our previous constraints and also consistent 
with each other. 

On the other hand, this approach is useful to verify the as- 
sumptions on the shape of the distribution functions used in 
Section l6T2l This is especially worthwhile for the ERDF, because 
the shape of the BHMF is relatively well determined at the high 
mass end, with the main uncertainty in the low mass slope, while 
the shape of the ERDF is poorly determined. Therefore, we as- 
sumed the double power law with fixed high mass slope parame- 
terisation of the BHMF. As shown above, the shape of the binned 
ERDF is consistent for all four assumed BHMFs, based on the 
4 best fit models. Thus it is justified to use one of these for the 
investigation of the ERDF shape. 

We again fix the break of the double power law and thus the 
only free parameter left is the low mass slope q-bh- We deter- 
mined the Eddington ratio selection function, using Eauation[T9l 
for a variety of values for q-bh, covering the whole range of ac- 
ceptable values. We use q'bh = -0.7 as lower limit, taken from 
the uncorrected BHMF, and q'bh - -3.2 as upper limit, corre- 
sponding to a single power law BHMF. The fitting results on 
these ERDFs with a Schechter function are given in Table 16.31 
While the normalisation of these ERDFs changes significantly 
for different assumed values of q'bh, the shape is not strongly af- 
fected and is consistent with our previous constraints thoughout 
the whole range. In particular, the ERDF is well described by a 
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log [Mq] log A 



Fig. 11. Same as Fig.[8]with the constraints from the 1 / Vmax method added (see Section l673] ). The binned results for the BHMF (left 
panel) were derived assuming the ERDF of the best fit solutions for the 4 models given in Table[3] Shown are only the two distinct 
models, the Schechter function (blue circles) and the log-normal distribution (green squares). Analogously, the binned results for 
the ERDF were derived assuming the BHMF of these 4 models. Shown are again only the two distinct models, the double power 
law (blue circles) and the modified Schechter function (green squares) 
. The binned results for the different models are consistent with each other, as well as with the results of Section [6^ shown as lines and by the 

shaded area. 



Table 4. Fitting results for the ERDF, determined using an ap- 
propriate Eddington ratio selection function, assuming different 
values for the low mass slope obh of the BHMF. 



Q-BH 



OA log A , x^/i-o-f- 



-0.7 


-1.51 


-0.63 


1.81 


-1.2 


-1.50 


-0.66 


1.77 


-1.7 


-1.51 


-0.69 


1.73 


-2.0 


-1.55 


-0.70 


1.73 


-2.2 


-1.59 


-0.71 


1.75 


-2.2 


-1.59 


-0.71 


1.75 


-2.7 


-1.79 


-0.73 


1.92 


-3.0 


-1.94 


-0.72 


2.02 


-3.2 


-2.04 


-0.71 


2.02 



Schechter function. While there is an indication for a flattening 
at the low A end, no indication for a real turnoff of the ERDF 
is present, as also shown in the right panel of Fig. [TT] A log- 
normal distribution is also appropriate, but needs to be cut off 
close to the maximum of the distribution. Thus it does not indi- 
cate a turnover, but only a flattening of the ERDF. We again want 
to emphasise that no prior assumptions on the ERDF are used 
here, we just modified the selection function using an assumed 
BHMF over a wide range of possible parameters. This strongly 
confirms our previous results for the shape of the ERDF, in that it 
shows that a Schechter function provides a good representation 
of the data. 



7. Discussion 

7.1. Active fraction of local black holes 

For a census of active black holes, the derived mass function 
of active black holes should be compared to the local mass 
function of quiescent black holes. Because the number of dy- 
namically measured black hole masses is still very low and 
the sample is inhomogeneous, the quiescent black hole mass 
function has to rely on the known M^h - bulge property rela- 
tions, thus converting galaxy luminosity or velocity functions 
into a black hole mass function. This approach h as been used 



by several authors to derive a local B HMF (e.g. Salucci et al 
'1999; ,Yu & Tremaine 2002; IShankar et al. 2004; Ma rconi et al 



20041). However, there is still some uncertainty in the estima- 
tion (iTundo et al.ll2007|). We com pare our active BHMF to the 
BHMF presented by iMarconi et al. (2004), shown as the solid 
line in Fig. [12] Our best fit model of the reconstructed active 
BHMF, derived above, is indicated as dashed line in Fig. [12] 

At this point we need to recall that our operational defini- 
tion of 'active' black holes only includes type 1 AGN. We are 
not able to distinguish between a true quiescent black hole and 
an AGN not selected due to obscuration. By dividing our active 
BHMF by the quiescent BHMF we thus get the fraction of black 
holes in an active stage, not hidden to our survey by obscuration, 
and thus a lower limit to the true active fraction. 

The lower panels in Fig. [T2]show the fraction of local black 
holes in an active stage as a function of the black hole mass, thus 
the black hole duty cycle. 

As circles we give the active fraction, or duty cycle, derived 
from the binned uncorrected BHMF, presented in Section 15.11 
where open symbols indicate bins that are affected by sample 
censorship. The estimate of the active fraction for the intrinsic 
BHMF is shown as dashed line, thus showing the intrinsic un- 
derlying black hole duty cycle. A clear decrease of the active 
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Fig. 12. Comparison of our active black hole mass function with the inactive BHMF of iMarconi et al. ] ( 120041) (soUd Une m upper 
panels). The circles show the uncorrected binned data, where open symbols correspond to bins that suffer from selection effects. 
In the lower panels the active fraction of black holes in the local universe is shown in logarithmic units. In the left panel the active 
BHMF and active fraction for the whole sample are shown. In the right panel the active BHMF and the corresponding active fraction 
are shown for two Eddington ratio bins (above and below log A = -1). The BHMF and active fraction for the best fit model of the 
intrinsic BHMF and ERDF are shown as dashed lines, in the left panel for the whole sample and in the right panel for the two A bin 
subsamples. There is a decrease in the active fraction with increasing black hole mass, in agreement with the cosmic downsizing 
picture. This decrease is already visible in the high A subsample (red squares). 



fraction with increasing Mbh is visible, being close to a power 
law with slope qaf ~ -0.86 over the whole covered mass range. 

Using a very different approach, IShankar et al.l (l2009l) pre- 
dicted the black hole duty cycle. They used simple black hole 
growth models, based on the local quiescent BHMF and the 
bolometric AGN luminosity function. They made the simpli- 
fied assumption of a single constant accretion rate, in contrast 
to the wide accretion rate distribution we assumed. Their active 
fraction also refers to the whole AGN population, while we are 
restricted to type 1 AGN. Nevertheless, when comparing their 
low z duty cycle with our results, we find an excellent agree- 
ment between both. However, taking into account th e large dif- 
ferences between the simple model of IShankar et al.l (|2009|) and 
our empirical determination, this agreement might even be a co- 
incidence. 

In the right panels of Fig. [12] we split our sample into two 
subsamples, based on the Eddington ratio A , at log A = -1. For 
both subsamples we computed the uncorrected BHMF and the 
active fraction. The uncorrected BHMF and the uncorrected ac- 
tive fraction for the low A subsample are shown as blue circles, 
while the blue dashed line shows the active fraction derived from 
the reconstructed BHMF (best fit modified Schechter function). 
Incompleteness sets in around 10^ Mq and is dominant below 
10^ Mq, therefore no information on the behaviour of the active 
fraction can be gained from these low A black holes. 

The high A subsample is shown as red squares, while the 
red dashed line shows the active fraction derived from the recon- 



structed BHMF, with the normalisation derived from the fraction 
of objects above log A = -I. The subsample is almost complete 
up to ~ 10^ Mq, where the low A subsample is already heav- 
ily incomplete. Above 10' the binned active fraction is in 
good agreement with the reconstructed intrinsic active fraction. 
This provides a consistency test for the reconstructed BHMF and 
ERDF estimate. But even without this comparison there is a clear 
trend present for the high A subsample with a decrease of the ac- 
tive fraction with increasing black hole mass, directly verifying 
our previous result from the uncorrected binned data. Thus, far 
more low mass black holes in the local universe are in an active 
state than high mass black holes. 



This result is in general agreement with the picture of 
anti -hierarchical growth of black holes (e.g. iMerlonil l2004t 
iMerloni & Heinzl 120081) . where the most massive black holes 
grew at early cosmic times and are preferentially in a less ac- 
tive stage in the present universe, and at present mainly smaller 
mass black holes grow at a significant rate, also known as cosmic 
downsizing. Our results strongly support this anti-hierarchical 
black hole growth scenario. This is in g eneral agreement with 
previous findings on low redshift AGN jHeckman et al] l2004t 
iGreene & H"oI|2007I: iGoulding et"ani20I0h that also report a de- 
crease of the active fraction for the most massive black holes, 
as we ll as with results at higher redshifts (IVestergaard & Osmerl 
I2009h . 
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Fig. 13. Comparison of the local BHMF of the HES with the 
BHMF presented in iGreene & Hoi ( l2009h . The blue asterisks 
and th e blue dotted line show th e BHM F from iGreene '&Hol 
(l2007h (corrected in IGreene & Hoi (l2009t) ). The open, black cir- 
cles show our BHMF, not corrected for evolution and sample 
censorship, while the filled, red circ les show our BHMF , with 
the black hole mass estimated as in IGreene & Hoi (l2007l) . The 
dashed line indicates our reconstructed BHMF for reference. 



7.2. The active black hole mass density 

We now want to estimate the black hole mass density of active 
type 1 AGN in the local universe 



Pact 



= I M.(^(M.)dM. , 



(21) 



with Myr 



using our results for the active BHMF. A 



lower limit on the local mass density of active black holes is 
given by the BHMF without a correction for sample censorship. 
We derived a lower limit of pact - 277 M© Mpc""'. Using our re- 
constructed BHMF, the local mass density of active black holes 
with log /I > -2 is then p^ct ~ nOOMoMpc"^^, a factor of 6 
higher then derived from the uncorrected active BHMF. The re- 
sults for the individual models are given in Table [3] 

The observational estimate of the integrated mass den- 
sity of the total black hole population i n the local universe 
is ptot ^ (3.2 - 5.4) X lO^Mp Mpc' ^ jShankar et al.1 120091: 
iGraham & Driveijl2007HYu & Lull2008l). U sing a value of 4.6 x 
10^ MoMpc"^, as presented bv lMarconi et al. (2004), results in 
a fraction of ~ 4 x 10"^ of the black hole mass that is currently 
actively accreting at a rate larger that 1 % of the Eddington limit 
(~ 6 X 10 for the uncorrected BHMF). 

7.3. Comparison with other surveys 

IGreene & Hoi (l2007l) presented a determination of the active 
black hole mass function for z < 0.352, using the SDSS DR4 
main galaxy sample as well as the QSO sample. They con- 
structed their sample based on spectroscopic confirmation of 
broad Ha lines, ending up with 8728 objects. For these they 
computed black hole masses from the Hff FWHM and line lu- 
minosity. 

As already mentioned in Paper I, an error has been discov- 
ered in the determination of the Vmax values in the work of 
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Fig. 14. Comparison of our local active BHMF (filled circles 
for uncorrected and dashed line for intrinsic BHMF) with the 
BHMF of the BQS, as determined directly as binned estimate by 
IVestergaard & Osmeil (l2009|)rblue asterisks) and as determined 
by a Bayesian approach bv lKellv et al.l (l2009l) (blue solid line for 
median and dashed lines for Icr uncertainty). 



IGreene & Hoi (l2007l) (J. Greene, private communication), result- 
ing into an erroneous luminosity f unction as well as BH MF. This 
error has recently been resolved (IGreene & Ho"2009^. Thus we 
cautio n not to use the original active BHMFjfrom Greene & Ho! 
(120071) . In Fig. [13] the active BHMF by IGreene & Hoi tQQ% 
from their SDSS sample is shown as blue asterisks. 

IGreene & Hoi (l2007l) have not taken into account the selec- 
tion effects caused by the use of the luminosity selection function 
and thus underestimate the number of active black holes at low 
masses, due to the discussed sample censorship. They also do 
not correct for evolution within their z ran ge. However, a direc t 
comparison with the mass function from IGreene & Hoi (I2009h 
can be made using our BHMF, without correction for evolution 
and sample censorship. 

For consistency, we also re-estimated the black hole masses 
of our sample, using the same formula as IGreene & Hoi (I2007h . 
using Ha FWHM and Ho- line luminosity. For our sample, the 
black hole mass distribution is shifted by 0.54 dex towards lower 
mass in the mean. Compared to the FWHM based Mbh this shift 
is 0.42 dex, thus ~ 0.1 dex can be attributed to the difference 
between the FWHM and o-\i„e based Mbh ■ The main reason for 
the remaining difference originates from a different relation of 
Ha luminosi ty to L^too found for our sample compared to the 
one given m IGreene & Hoi (l200l . This difference leads to an 
offset of 0.31 dex. The remaining offset can be attributed to the 
different Rblr - L scaling relation as well as to scatter in the 
relation between the FWHMs, as shown in Fig.[T] 

The resulting BHMF of the HES is shown as filled circles 
in Fig. [13] Both BHMFs are fully consistent with each other, 
especially at the high mass end, where different survey selec- 
tion effects are not important. At the low mass end the SDSS 
BHMF seems to exhibit similar survey selection effects as our 
HES sample, resulting in a consistent uncorrected BHMF, even 
at the biased low mass end. 

Recentlv. Yesterg aard & Osmerl (l2009l) presented the binned 
local active BHMF o f the Bright Quasar Survey (BQS, 
ISchmidt&Greenlll983h . in the redshift interval z - - 0.5. In 
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Fig.[T4]we compare their derived BHMF with our binned BHMF, 
not corrected for evolution and sample censorship (filled black 
circles), and our reconstructed intrinsic BHMF (dashed black 
line). We also show the local BHMF of the BQS a s blue soUd 
line, but determined using a Bayesian approach (iKellv et alj 
l2009h . 

The most direct comparison between the BQS and the HES 
is with the binned estimates. At the high mass end, both binned 
estimates are in reasonable agreement. However, the BQS does 
not cover exactly the same redshift range as our HES sample. 
This might also cause some difference between both BHMFs, 
due to evolution of the BHMF, which has the largest effect 
at the high mass end. Because the BQS is not as deep as the 
HES, incompleteness sets in at higher Mbh in the binned BHMF. 
Also it is known that the B QS suffers fr om bot h incompleteness 
dGoldschmidt et al.l [19921 iKo ehler e t al.||1997h as well as over- 
completeness dWisotzki et ani2000l) . Thus, the HES is superior 
to the BQS f or a determination of the local active BHMF. 

Recentlv. lKellv et al] (l2009h presented a determination of the 
active BHMF from the BQS using a Bayesian method, taking 
also into account scatter in Mbh and accounting for black holes 
below the flux limit of the survey. Their approach aims at cor- 
recting their BHMF for sample selection effects caused by the 
flux-limit, as we did in Section |6l However, they modeled the 
BHMF with a combination of Gaussian functions and also en- 
forced a log normal distribution for the ERDF, while we mainly 
used a Schechter function description without a specific maxi- 
mum and with a high fraction of objects at low A . In Fig. [14] we 
compare their posterior median BHMF (blue solid line) with our 
intrinsic BHMF (black dashed line). While both mass functions 
are consistent at the high mass end, there is a clear disagreement 
at the low mass end. Their BHMF is rather consistent with our 
uncorrected BHMF. We speculate that the reason for this dis- 
agreement might lie in the different assumptions on the shape of 
the BHMF and ERDF. This emphasises the importance of the as- 
sumed ERDF for the determination of the underlying BHMF. An 
important constraint on the ERDF is provided by the condition 
to recover the observed luminosity function as a convolution of 
BHMF and ERDF, as we have ensured by construction. 

So far, little observational results exist on the distribution 



function of Eddington ratios from AGN surveysJ Vu et a l; ( 2005h 
used a sample of type 2 AG N from the SDSS jKauffma nn et alJ 
l2003t lHeckmanet"ani2004 to determine the ERDF T heir re- 
sults h ave recently been compiled by Hopkins & H ernquisll 
(120091) . They also fitted the ERDF by a Schechter function and 
found an average slope of oa ~ -1 .6. Our constraints on the lo- 
cal type 1 ERDF presented here are consistent with this average 
slope of the ERDF of type 2 AGN. This might indicate a similar 
accretion behaviour of type 1 and ty pe 2 AGN, as exp ected from 
the standard unification model (e.g. I Antonucci|[T993h . 



7.4. Evolution of the active fraction 



IVestergaard et al.l (l2008h presented a determination of the 
BHMF in the redshift range 0.3 < z < 5. They used a well- 
defined, homogeneous sample of 15180 quasars from the SDSS 
DR3, already used by Richards et al. (2006b) for the construc- 
tion of the luminosity function. They found a high-mass decline 
with constant slope /3 x -3.3 at all epochs. Our high-mass slope 
of j6 = -3.1 for z < 0.3 (when not corrected for evolution) is 
consistent with their higher-z result within the uncertainties. 

We compare our z < 0.3 BHMF, not corrected for evolu- 
tion within the z-bin and for sample censors hip, with the low- 
est redshift bin (0.3 < z < 0.68) BHMF of IVestergaard et alJ 
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Fig. 15. Comparison of our uncorrected z < 0.3 BHMF (filled 
circles) with the active BHMF within 0.3 < z < 0.7 from the 
SDSS ( IVestergaard et al.|[2008l) . shown as the blue dashed line 
and aster isks, and the active BHMF o f the LBQS within 0.2 < 
z < 0.5 (^Ve stergaard & Osmej|2009l) . shown as the red dotted 
line and triangles. 



(120081) . shown in Fig. [15] We also show t he active BHMF of the 
Large Bright Quasar Survey (L BQS; e.g. Hewett et al. 2001) in 
the redshift bin z = 0.2 - 0.5 (IVestergaard & Osmer 2009.) as 
triangles. Both SDSS and LBQS BHMFs are in general agree- 
ment, even though they do not cover exactly the same redshift 
range. The decline of the space density at the lowest Mbh in both 
BHMFs is mainly due to incompleteness in this mass range in 
the SDSS QSO sample as well as in the LBQS QSO sample. At 
the high mass end the BHMF shows a similar slope but a larger 
space density than our HES BHMF. This seems to indicate evo- 
lution of the BHMF between these redshift bins. Because the 
mass function of the total supermassive black hole population 
will only decrease at the high mass end toward higher z, we can 
use the local quiescent BHMF as an upper limit for the mass 
function at 0.3 < z < 0.68. This then implies an increase of the 
active fraction at the high mass end towards higher redshift, ex- 
actly as would be expected in the cosmic downsizing scenario. 
Thus, the number of the most massive black holes being in an 
active stage in the present universe seems to be lower than at 
earlier cosmic epochs. 

8. Conclusions 

We have presented a study of the low-redshift active black hole 
population, residing in broad-line active galactic nuclei. We es- 
timated black hole masses and Eddington ratios, and from it 
estimated the local active black hole mass function and the 
Eddington ratio distribution function. Our sample was drawn 
from the Hamburg/ESO Survey and contains 329 quasars and 
Seyfert 1 galaxies with z < 0.3, selected from surveying almost 
7000 deg^ in the southern sky. 

We estimated black hole masses from single-epoch spectra, 
measuring the line dispersion of the broad HfS line and the con- 
tinuum luminosity at 5 100 A L5100, using the common virial 
method. We took care to avoid contamination of the line mea- 
surement by neighbouring emission lines and roughly estimated 
the degree of host galaxy contribution to Lsioo- This has been 
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found to be negligible for the most luminous AGN and not dom- 
inant even at the low luminosity end of our sample. We applied a 
rough statistical correction to the continuum luminosities to take 
into account the host contribution. The bolometric luminosity 
and thus the Eddington ratio A , has been estimated from Lgioo- 

The observed black hole masses cover a range 10^-2-10'' Mo 
and the Eddington ratio is roughly confined between 0.01 - 1. 
The observed distributions of these quantities are understood 
by the underlying distribution functions of black hole mass and 
Eddington ratio, in combination with the survey selection func- 
tion, as we explicitly demonstrated by Monte Carlo simulations. 

We made an attempt to determine these two distribution 
functions in an unbiased way. First of all, when we want to de- 
termine the active BHMF, we have to make clear what we mean 
by an active black hole, due to the wide distribution of accretion 
rates. We used a lower Eddington ratio cut of log A - -2, in 
agreement with the observed range of Eddington ratios. Using 
a different cut for A will preserve the shape of the BHMF, but 
change their normalisation, due to our assumption of an uncorre- 
cted BHMF and ERDF This is already shown in the left panel of 
Fig[T2] The normalisation and therefore the space density clearly 
depend of the chosen definition of an active black hole. 

Next we have to be aware that our sample is selected on AGN 
luminosity, not on black hole mass. Therefore, we have to make 
sure that we properly account for active black holes below the 
flux limit of the survey. We presented a method that determines 
the active BHMF as well as the ERDF at the same time, by a 
maximum likelihood fit. Here, the bivariate probability distribu- 
tion of black hole mass and Eddington ratio is fitted to the ob- 
servations. This probability distribution is given by an assumed 
BHMF, ERDF and the selection function of the survey. We also 
corrected for evolution within our redshift range, transforming 
the distribution functions to z = 0. This maximum likelihood 
method also ensures the consistency of the derived BHMF and 
ERDF with the AGN luminosity function. We were able to put 
tight constraints on both the active black hole mass function and 
the Eddington ratio distribution function. 

The Eddington ratio distribution function is well described 
by a Schechter function with low A slope ff^ ^ -1.9. The data 
are consistent with no decrease of the ERDF at low A , within the 
constrained range. Using a log-normal distribution, we found a 
maximum at log /I = -1.8, what can be taken as an upper limit 
for a potential turnover in the ERDF. Our results clearly show 
a wide distribution of Eddington ratios, in contrast to a single 
value or to a narrow log-normal distribution, which is based on 
the observed distribution, without accounting for the underlying 
selection effects. While we also observe a narrow log-normal 
distribution of Eddington ratios, this is in agreement with the 
constrained Schechter function or wide log-normal distribution 
for the Eddington ratio distribution function, when survey selec- 
tion effects are properly accounted for, because \ow-A objects 
will be systematically missed in flux limited samples. 

The active BHMF is well described by different analytic 
models. In general, it strongly decreases at the high mass end 
and follows a power law at the low mass end with slope of 
a ^ -2. A good fit to the data is achieved by a function similar 
to a Schechter function, but modified by an extra parameter that 
determines the steepness of the high mass decrease. We found 
no evidenc e for a decrease of th e BHMF toward low mass, as in- 
dicated by [GreenT&H^ dZOO^) for Mbh < IO^-^Mq. However, 
our sample is not very sensitive in this low mass range. 

We compared our local active BHM F with the local quies- 
cent BHMF from iMarconi et al] ( |2004|) . determining the active 
fraction, or duty cycle, of local black holes. This active frac- 



tion is decreasing with increasing black hole mass, consistent 
with a power law with slope ~ -0.86. Thus, the most massive 
black holes in the present universe are less active than their 
lower mass companions. At the highest MeHonly 10""* of all 
black holes are currently in an active stage, i.e. accreting above 
0.01 of the Eddington rate. This supports the general picture of 
anti-hierarchical growth of black holes. This mass dependence 
of the active fraction indicates that our assumption of an uncor- 
rected BHMF and ERDF cannot be sustained up to low values 
of A and thus we caution to extrapolate the distribution functions 
into the low A regime. Investigating a mass dependence of the 
ERDF would especially require a wider luminosity coverage of 
the sample. 

By comparing our low z BHMF with the BHMF of 
a highe r z-bin, p resented by IVestergaard et al.l (l2008h and 
Vesterga ard & Osm er (2009|), we found an indication that the 
most massive black holes are currently in a less active stage 
than at earlier cosmic times, also in general agreement with anti- 
hierarchical black hole growth. 

Recentlv. lMarconi et al.l (l2008h proposed a modified method 
to estimate Mbh , taking into account the effect of radiation pres- 
sure. So far, it is still unknown if radiati on pressure h as an impor- 
tant effect on the BLR or not (see e.g. lNetze3l2009l) . If we take 
into account radiation pressure and apply their Mbh estimation 
formula to our sample, the major effect is an increase of 
Mbh especially for the low Mbh objects. In total, the dispersion 
of the Mbh distribution decreases from 0.65 dex to 0.63 dex. In 
the BHMF the space density at median Mbh increases, while 
at high Mbh the space density slightly decreases. This would 
strengthen even further the evidence for anti-hierarchical black 
hole growth. On the other hand it would change our observed 
Mbh , and especially our A , distributions, and thereby our con- 
strained BHMF and the Eddington ratio distribution function. 

Our work strengthens the scenario of a nti-hi erarchical 
growth of black holes , also seen in other studies (iMe rloni 2004 ; 
Heckman et al 1l200l iGreene & Holl2007t IShankar et al.. .20091: 



Vestergaard & Osmeijl2009l) . at least at low redshift. The obser 



vation of 'cosmic downsizing' in the X-ray luminosity function 
(e.g. lUeda et al . 2003: Hasinger et al. 20()5i), as well as in the 
optical, radio and IR luminosity function (e.g. Hunt et al] 20041 : 
ICirasuolo et aP 120051: iMatute et al.l l2006t ICroom et al.1 l2009l) . 
i.e. the flattening of the faint end slope of the luminosity function 
towards higher redshift, is explained by the shift of the typical 
black hole mass of an active accreting black hole toward lower 
mass. 

The presented local active black hole mass function and 
Eddington ratio distribution function serve as a local anchor for 
future studies of both distribution functions. These will provide 
further information on the cosmic history of growth and activity 
of supermassive black holes. 
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